Critical Temperature of Dilute Bose Gases 



Volker Betz and Daniel Ueltschi 
Department of Mathematics, University of Warwick, Coventry, CV4 7AL, United Kingdon^\ 

We compute the critical temperature of Bose-Einstein condensation in dilute three-dimensional 
homogeneous Bose gases. Our method involves the models of spatial permutations and it should be 
exact to lowest order in the scattering length of the particle interactions. We find that the change in 
the critical temperature, AT C /Tc°\ behaves as cap 1 / 3 with c = —2.33; this contradicts the current 
consensus among physicists. 

PACS numbers: 05.70.Fh, 03.75.Hh, 05.30.Jp 
Keywords: Dilute Bose gas, Bose-Einstein condensation 



1. Introduction. The study of the effects of parti- 
cle interactions on the critical temperature of the Bosc- 
Einstcin condensation has a long and tortuous history. 
Several studies found that repulsive interactions decrease 
the critical temperature [Iml, but most studies pointed 
to an increase of the critical temperature [4H12I]. The 
current consensus among physicists is that the change 



(a) 



Tc°^ behaves as 



r. 



(0) 



cap 



1/3 
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with positive constant c ps 1.3 [9H12J. Here, de- 
notes the critical temperature of a dilute homogeneous 
Bose gas and a is the scattering length of the (repulsive) 
potential that describes the par ticle interactions. These 
results are reviewed in Refs [13l - ll5l |. where more details 
and additional references can be found. The latter article 
[lH proposes a mathematically rigorous upper bound on 
the interacting critical density; it is useful but it does not 
settle the question. 

The present article uses a different approach that in- 
volves "spatial random permutations". We initially ex- 
pected that our calculations would confirm the consensus 
among physicists, but they rather seem to invalidate it: 
We find a negative constant c = — 3^(3/2)4/3 ~ —2.33, 
which implies that interactions decrease the critical tem- 
perature. While not rigorous, we expect our approach to 
be exact and to yield the correct constant. But we also 
trust the physics literature, and we are left puzzled. We 
do not know how to resolve these contradictions. 

2. Overview of the method. We start with 
the Feynman-Kac representation of the Bose gas, where 
quantum particles become winding Brownian bridges in 
one more dimension 0, [l6j]. For dilute gases, Bose- 
Einstein condensation is signaled by the occurrence of 
infinite loops [13, EH ■ We actually consider the simpler 
models of spatial random permutations where Brownian 
bridges have been integrated . We derive a simplified 



model where the original two-jump interactions are re- 
placed by "cycle weights" . We argue that the new model 
has almost the same marginal distribution on permuta- 
tions, and the same critical temperature to lowest or- 
der in the scattering length. The critical temperature of 
the model with cycle weights can be computed explicitly. 
Many of the concepts and methods in this article can be 
fully justified mathematically. There are five instances, 
however, where a leap of faith is required. We have out- 
lined these steps by writing them as "conjectures" . While 
the validity of the conjectures 1, 2, 4, and 5, seems clear, 
the conjecture 3 is tricky. If our calculations do not give 
the correct result, the most likely reason is the failure of 
this conjecture. 

Finally, we put our method to a test by computing the 
free energy of the effective model. As it turns out, it is 
equal to the free energy of the dilute Bose gas [2(J, Hl[ , 
to leading order in a. This does not prove the method 
right, but this raises our confidence in it. 

3. Mathematical setting. The state space for N 
bosons in a box A C I 3 is the Hilbert space Lf ym (A N ) 
of square-integrable complex symmetric functions. The 
Hamiltonian is the Schrodinger operator 
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with Ai the Laplacian for the i-th particle, and U a func- 
tion that acts as a multiplication operator. We always 
suppose that U is nonnegative, spherically symmetric, 
with finite range. It is characterized by its scattering 
length a. If U is small and integrable, a « J U(x)dx 
is the first Born approximation to the scattering length. 
If U consists of a hard core, a is the radius of the hard 
core. We refer to Appendix C of the monograph [22| for 
the general definition and more information. 

4. Feynman-Kac representation. The Feynman- 
Kac formula allows to represent quantum particles by 
Brownian bridges [3, EH . The partition function of the 
Bose gas can be written as that of a model of "spatial 
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permutations" [lj|, namely, 
z = Tr i2 ym (A«) e _/3ff 

m(4*p)w/* ^ J AN axe 

where Sn denotes the set of permutations of N elements, 
x = (xi, . . . , xn) G A N , and where the "Gibbs factor" 

C -H 1 (x : ir) j g gj ven 

N . 
i— 1 

r 1 / >2/3 -i 
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In this article, W* „ denotes the Wiener measure for the 
Brownian bridge traveling from x to y in time t. It sat- 
isfies 

jdW t Xty (u } )=g t (x-y), (5) 

with gt the Gaussian function, namely, 

g t (x) = {2Kt)-V 2 e-* 2 l 2t . (6) 

We temporarily ignore boundary effects, but we will need 
to take them into account later when we calculate the 
weights of macroscopic cycles. 

Let Ho(x,ir) denote the Hamiltonian of the model of 
spatial permutations that corresponds to the ideal Bosc 
gas: 

1 N 

h o(x,tt) = —^{xi-x^l 2 . (7) 

The equation ((4]) for Hi can then be written as 

N r 

»=1 •* 

r 1 r 2P i 

cx p{-2 ^ / U ("i( s )-"i( s )) ds }- (8) 

Here, we introduced the normalized Wiener measure 
W*^, that is equal to g t _1 (x - y)W* y . 

5. Cycle lengths and Bose-Einstein condensa- 
tion. The order parameter for Bose-Einstein condensa- 
tion is Penrose and Onsager's "off-diagonal long-range or- 
der" [23[ ■ But following Siito [ijj , we rather consider the 
length of permutation cycles. It is expected that Bose- 
Einstein condensation is accompanied by the occurrence 
of infinite cycles. 

CONJECTURE 1. The critical temperature for Bose- 
Einstein condensation is identical to the critical temper- 
ature for the occurrence of infinite cycles. 



This conjecture is compatible with the conclusions 
drawn in Ref. [18]; it was indeed argued that it is valid 
for dilute systems, but that it may be invalid otherwise. 
From now on, we rover in the realm of spatial permuta- 
tions, and we invoke methods and heuristics of classical 
statistical mechanics and of probability theory. The first 
step consists in deriving a model with "two-jump inter- 
actions" . 

6. Model with 2- jump interactions. We consider 
the Hamiltonian 

H 2 (x,ic) =H (x,w)+ v ij( x >n) (9) 

where Vij{x, it) depends on Xi, x^u), Xj, x n (j), and it rep- 
resents the interactions between the jumps Xi H ► x-xU) 
and Xj i — y x^fjy It involves the particle interactions and 
it is given by the formula 

( 10 ) 

The general properties of Vy are not immediately appar- 
ent. By the Feynman-Kac formula, it is related to the 
integral kernel of the operator e 2/3A — e^ 2A ~ c/ - ) . Let us 
notice that, when U is a hard core potential of radius 
a, Vij is equal to the probability that a Brownian bridge 
from Xi — Xj to x v u\ — 2^r(j) intersects the ball of radius 
a around the origin. 

The formula ([TP]) was computed in Ref. [24], where it 
was also shown that, in a certain sense, 

ffi(x,7r) = H 2 (x,tt) + 0(a 2 ). (11) 

An expression for 0(a 2 ) can be found in [24[ that results 
from a cluster expansion. Rigorous estimates for this 
term have yet to be derived, though. 

CONJECTURE 2. The critical temperature of the 
model H% is identical to that of the model H\, up to a 
correction o(a). 

While a bit simpler, the new model is still intractable 
per se and further simplifications are necessary. The 
next goal is a model of spatial permutations with "cy- 
cle weights" , see Eq. (l23l) below. 

The correction 0(a 2 ) in Eq. (fTTj) also involves the in- 
verse temperature fj and the density p. One should dis- 
cuss the dimensions that are present here. H\ and H 2 
are dimensionless; a is like a length; fj is like the square 
of a length; p^ 1 is like the cube of a length. It seems 
that the correction in pip is really like 0(((jpa) 2 ), which 
is dimensionless. It is not clear at this stage whether it 
is enough to suppose that fj pa 1, or whether an ad- 
ditional condition is needed. In this article we always 
suppose that fjp 2 l z ~ 1, i.e. the system is in the regime 
of the phase transition. We only characterize the error 
terms O(-) and o(-) by their dependence on the scattering 
length a. 

7. Using spatial averaging. We now describe the 
key step that allows to replace the 2-jump interactions 
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Vij(x,w) by cycle weights ai within cycles of length £, 
and aijt between cycles of length £ and £'. Cycle weights 
are much simpler because they do not depend on spatial 
positions. The difficulty is to make this simplification 
while keeping the critical density unchanged, at least to 
lowest order. 

Let 9 be a random variable that depends only on the 
permutation, not on the positions. Its expectation can 
be written as 



E{6) = i- J dxJ2^ 



(12) 

We introduced the probability measure fj,^' that depends 
on 7r, and which is defined by 



d^>(x) 



1 



(13) 



with 



Z w = / e' Ho{x ^dx. (14) 

We now introduce a "permutation free energy" by 

e- F W = /" d^M^e-^ ^ 3 ^ . (15) 

The main step of our method is to take the integral inside 
the exponential, namely 

c- F M w exp{ - / d^ W (a;) V U (*>*")}• (16) 

The left side is larger than the right side by Jensen's 
inequality, but it is hard to formulate a converse inequal- 
ity. Let us carefully justify Eq. (fll)|). We need to take 
the large size of the system into account, and to think in 
terms of typical positions. For a given permutation, let 
y be a typical realization of the measure fi^\ and let z 
be a typical realization of the measure 
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-H (x,tt) e -J2i <j Vij(x,Tr) 



normalization 
For all macroscopic observables A we expect that 

A(y)=A(z)(l + 0(a)). (17) 

This holds in particular when the macroscopic observable 
isY*<jVij(x,ir). We have e~ F ^ « e "^«<i y « (z,7r) . 
We replace it by e~ F( -^ « e~ v «(f' ,r ) 1 the difference 
of free energies should be of order 0(a 2 ). The latter gives 
the expression p^|) . 

We get the new Hamiltonian 

H 3 (x,w)=Ho(x,k) + Y^ f V l3 (y,7T)d^\y). (18) 



CONJECTURE 3. The critical temperature of the 
model i?3 is identical to that of the model H2, up to a 
correction 0(0). 

8. Model with cycle weights. The expression 
(fT8|) can be simplified further. If i and j belong to the 
same cycle 7 of length i, we have 



V ij (y,n)d fJ ,M(y) = I V lJ (y, 1 )d^>(y), (19) 



where the new expectation refers to a system of £ points 
with the cyclic permutation 7. If j and j belong to dif- 
ferent cycles 7 and 7' of respective length £ and £', we 
have 



V ij (y,n)d^\y)= / V i3 (y, 7 U 7') d^ (7U7,) (y), 



( 2 °) 

where 7 U 7' denotes the permutation of £ + £' elements, 
with two cycles of length £ and £'. Let us introduce 



ae,e< - , 



J2 Jv iJ (y, 1 )d^\y), (21) 

\ E /^(y,7U7')dM (7uy) (y)- (22) 



Notice that these expressions depend on the cycle 
lengths, but not on the explicit cycles. We get a Hamil- 
tonian for spatial permutations with "cycle weights" , 
namely 



H 3 (x,ir) =H (x,tt) +^2(u£- ai t i)ri(n) 



> 1 



(23) 



> 1 



An important remark is that the weights an and <xi t gj 
depend on the domain A of the system. For the small 
cycles, this dependence is anecdotal as it disappears in 
the infinite volume limit. For macroscopic cycles, con- 
finement in a bounded domain is important; since the 
cycles cover the domain, all points find themselves at fi- 
nite (spatial) distance of other points that are very far 
away along the cycle. We need not worry about cycles of 
intermediate size, since these are known to involve only 
a vanishing fraction of particles Qj], Oil HH . 

9. Simplified 2- jump interaction. To first order, 
the 2-jump interaction is given by 



1 



4/3 



dw: 



4rf 



V(xi Xj,X^(j^ 3'7r(j))' 



(24) 



The notation = means that the equality holds to first or- 
der in a. Recall that the first Born approximation to the 
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scattering length is g^- J U(x)dx = a. By the property of 
the Wiener measure, 



dW^(w)U(u;(s)) = / dzU(z) 



g s (x - z)g i0 - s (y - z) 



34/3 - y) 



It follows that 



(25) 



V(x,y)=2na r ^^-^ ds. (26) 

Jo 94/3{x-y) 

The following formula about Gaussian functions is not 
hard to establish, and it is useful here: 

r— = t g st (t- s ) {{t - s)x + sy) . (27 

g t {x-y) 

Then, after a change of variables, 

V(x,y) = 8n(3a / 3 4 /3s(i-*) ((1 - s)x + sy)ds. (28) 
Jo 

10. Weights of finite cycles. We now calculate 
the one-cycle weights to first order in a. We do it first 
for finite cycles, which will lead to the result in Eq. (131)1) . 
and we will do it in the next paragraph for macroscopic 
cycles. Recall the definition (f2"Tj) of the weight on. We 
have 



/ ^-(y,7)dM (7) (2/)- 

3=2 J 



(29) 



Gaussian functions satisfy g s * gt — g s +t- We use this 
equation repeatedly, in conjunction with (fTB")) and (f2"5|). 
and we get 

1 r 1 

a e = (47r/3£) 5 / 2 aV / F^(s)ds, (30) 

,-9 JO 



with 



F tj( s ) = TT7 / dxid^dxjdxj+i 32/3(^1 - ^2) 

l A l JA 4 

92p{j-2)(x2 - Xj)g 2 p(xj - a^+Offa/a^-^fe+i - £i) 
34/38(1-8) ((1 - s)(x 1 - Xj) + s(x 2 - x, J+1 )). (31) 

We used the fact that = \A\g 2 pt(0) = |A|(47r/3fT 3/2 - 
We have neglected boundary effects, but they are irrel- 
evant for finite I. We now consider the case j = 2 and 
I ^ 3. We actually have 

Fn{s) = 777/ da; 1 dx 2 da;3 32/3(^1 - x 2 ) 
l A l JA 3 

92p{x2 - x 3 )g 2 p(i- 2 )(x3 - Xi) 
34/3s(i-8) ((1 - s){xi - x 2 ) + s(x 2 - xi)). (32) 



Using translation invariance, we can set x 2 = 0, picking 
up a factor |A|. We get 



F t2 (s) = (4nP)- 6 (£ - 2)- 3 / 2 (2s(l - s)) 



-3/2 



da; 1 da; 3 



exp 



(x 3 -xi) 2 ((1 - s)xi - sx 3 ) 2 



Af3 4/3 A/3(£ - 2) 



8/3s(l - s) 



(33) 



The space variables factorize with respect to space direc- 
tions. Integrating the Gaussians, we get 



Ft 



I2\ 



(s) = (4^)-V - 2)- 3 / 2 (2 S (l - s) )- 3 / 2 (-g=) 

(34) 

with the matrix A given by 



A = 



1 + — 4- — 

1 ' 1-2 > 2s 

_1 1_ 

2 t-2 



_1 1_ 

2 £-2 

1 + 7=2 + 2(l-s) 



(35) 



The determinant of this matrix is equal to (2s(l — 
•■OrH 1 + 1=2)- We obtain F t2 {s) = (An/3)- 3 (£ - I)- 3 / 2 . 
Let us note that it does not depend on s. This result 
holds for I ^ 3 and for j = 2, and also for j — £ which 
is the same. We now treat the case £ ^ 3 and 2 < j < £. 
We introduce z\ = x\ — Xj and Z2 = X2 — Xj+i in (|3 1 : 
we replace x 2 by using translation invariance, and we 
set X\ — x. We obtain 



F tj (s) = (An0)-^((j-2)(e-j)) 6/2 (2s(l-s)) 6 < 

f c x 2 (zt- x) 2 (x-z 1 + z 2 ) 2 
J 4$ - ~p 



(x + z 2 ) 2 _ ((1 - a)zi - sz 2 f 
i/3{£-j) 



3s(l - s) 



} 



= (4^)-/ 2 [( J -2)(^)(2 S (l- S ))]- 3/2 (^) 3 - 

(36) 

Here, the matrix A is given by (with u = 1 + an d 
«=1 + A) 




(37) 

One easily finds detA = uv/2s(l — s). It follows that 
F tj (s) = (4tt/3)- 3 (j - ir 3 / 2 (^ - j + I)- 3 / 2 . (38) 

The formula for I 3 and j = 2,^ turns out to be 
compatible with the formula above. We have checked 
that it also holds in the case £ = j = 2. Notice that (s) 
does not depend on s. After minor rearrangements, we 
obtain 



£a s-^/ 



(4^)1/2 j^\j{£-j) 



3/2 



(39) 
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It is not hard to see that 
2(( 



lim 



Q'f 

T 



(4tt/3) 1 /2 



(40) 



This represents the energy per particle in a large finite 
cycle. 

11. Weights of macroscopic cycles. When the 
cycles are macroscopic, i.e. £ ~ N, we can no longer 
assume all points to be far away from the boundary of 
A; boundary effects have to be taken into account. For 
the Bose gas with periodic boundary conditions, (jB]) has 
to be replaced with the periodized version 



fit 



(A) 



(x) = 9t{x + Ly). 



(41) 



Here, L denotes the size of the box. The identity 
g% * = g\+ a still holds, where the convolution is 
now defined by / * g(x) = J A f(x — y)g(y) dy. Thus ([51]) 
holds with g replaced by g\ throughout, and indeed this 
is the correct finite volume expression even for finite cy- 
cles. For £ finite and A — > oo, we have g^eix) — > g2/3i(x)> 
which justifies our using g instead of g^ in d3U). But 
for £ — eN with any e > 0, it is not hard to see that 
l^lls^k 2 -) — > 1 for all x. Those Gaussians with macro- 
scopic variance in (|3"Tj) must be replaced by |A| _1 . A 
similar calculation as above then leads to the result 



a e = 8n(3p^a£ + - Ca/T 1/2 

|A| 



(42) 



The constant is equal to C = (3 — §7i/2)/-\/^> see Eq. 
()56|) . but this term does not play any role. 

12. Two-cycle weights. Recall the definition (|22|) 
of the weight a^y . We have also made computations in 
this case, and we have found that 



OtlJi 



4irf3££'a 



(43) 



In retrospect, we could have guessed this formula: From 
the definition (|2"2"1) . we clearly have 



a e ,t> = \££' 



\ J+l {y n yji)d^'\y). (44) 

Retracing our steps back to the Feynman-Kac representa- 
tion of the quantum model, we see that J Vi^+id// 71 " 17 ' 
gives the average interaction between the first Brownian 
bridge of the first loop and the first Brownian bridge of 
the second loop. Let us fix the first loop, and let us fix 
the imaginary time at which the interactions are com- 
puted. The integral over the second loop gives 87ra/|A|, 
after dividing by the normalization. This quantity does 
not depend on the first loop nor on the imaginary time. 
We get a factor j3 with the time integration, but the inte- 
gral over the first loop is canceled by the normalization. 
We then get fgg]). 



Because a^i is proportional to ££', there is an impor- 
tant simplification: 



E 

jt> i 



|A| 



(45) 



so that this term is a constant in the canonical ensemble. 
It does not affect the probability distribution of permu- 
tations, hence the critical density. We keep it, however, 
in order to perform the test of the free energy. 

Let us subtract the limit (|40j) in the one-cycle weights, 
and define 



2£a 
N) 



e-i 



We now consider the Hamiltonian 



3/2 



C(l) 



(46) 



|A| 



8TT0pi o) aN 
H (x,n) + ^2 v-itM- 



(47) 



H4 is essentially identical to H3, we only replaced certain 
expressions by their lowest order approximations. We 
should emphasize, however, that the approximating of 
aij: has substituted the interactions between different 
cycles by a constant, which is a huge simplification. 

CONJECTURE 4- The critical temperature of the 
model Hi is identical to that of the model H3, up to a 
correction o(a). 

13. Properties of one-cycle weights. We now 
show that a' t = Ca(l + O^" 1 / 5 )), see Eq. (56]) below. 
We will actually not use this formula directly, but it helps 
gaining a better understanding of the method. 

Let us simplify the notation a bit and introduce <j>i 
such that = 2(47r/3)~ 1/2 a<^. We have 



j(£-j) 



3/2 



1 

-372 



3>t/2 ■ 

We assumed here that £ is odd. For even £, we need to 
subtract a term 2 3 / 2 /\/^, which is smaller than the error 
term. The last term in (gSJ is equal to 2 3 / 2 V£+0{£- 1/2 ). 
For the first term, let us first observe that 



e/2 



lim 



I 



j(t-3) 



3/2 



1/2 



1 

~3/2 



(1 



\3/2 



- 1 



dx 



(49) 



Let £ > 0. We have 



i/2 

E 



[( 



= V~£ 



3(1-3) 
1/2 r 



3/2 



(1 - x y/ 2 



,•3/2 



- 1 



dx 
Z3J2 



+ 0{{e£y 1/2 ). (50) 
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By making the change of variables x = sin 9, and then 
using trigonometric identities, we find 



1/2 



d.r 



2-4e 



C 3/2 (1 _ x) 3/2 ^{T^T) 



(51) 



One also has J 1/2 x~ 3 / 2 dx = 2/^/e-2 3 / 2 . It follows that 



1/2 



1 



1 



dx 



.{1-x) 3 / 2 
Next, it is not hard to check that 

el 



2 3 / 2 -3Vi+0(£ 3/2 ). (52) 



3/2 



1 

372 



E 



0(e 2 £) 



= 3Vrf - 3 + § 7l/2 + 0(e 2 £) + 0{{e£)~ 1/2 ), (53) 

where 71/2 is a generalized Euler constant, that is defined 
by 



7l/2 



= lim 



r " 1 r r 



da; 



= 0.5396. 



(54) 



Collecting all the necessary terms, and choosing e — 
£- 3 / 5 , we obtain 



^ = -3 + f7i/2 + 0(^ 1/5 ). (55) 

It follows that 

a' e = -(6 - 3 7 i/ 2 )(4^)- 1 / 2 a(l + 0{£- 1 ' 5 )). (56) 

We observe that the weights a' e are negative. This is 
important, as it will lead to the decrease of the critical 
temperature. Notice that the sign of a'g was not obvious 
prior to the computation. 

14. The critical density with cycle weights. 
Let H§ be like H4 but without the constant terms. Ex- 
plicitly, 



(57) 



These Hamiltonians generate the same Gibbs distribu- 
tions, so that they have the same critical density. We 
now calculate the critical density of the model H$ using 
the grand-canonical ensemble. Namely, we consider the 
pressure 



FaT^E 



N > 



„ (47r/3)TriV! 



(58) 

It can be computed by going to the Fourier space and by 
introducing occupation numbers, as for the usual ideal 
Bose gas. In the limit of infinite volumes, we find 



p(/M= E 



< y (4w) 3 / 2 (P£) 5 / 2 ' 



(59) 



The pressure is finite for /1 ^ and it is increasing in 
/i. It is infinite for /1 > 0. As in the ideal gas, the critical 
density is given by the derivative of p at \i = 0: 



M=0- 



E 



' (4tt/^) 3 / 2 



(60) 



We need to emphasize that this formula for the criti- 
cal density is computed by analogy with the ideal Bose 
gas and its relevance for macroscopic cycles is far from 
clear. In the case oij = (where p c = pi " 1 ), Siito proved 
that it is indeed the critical density for the occurrence of 
infinite cycles [17j. This result was extended to the case 
of vanishing cycle weights — such that a' e — > faster 
than 1/log^ as £ — > 00 in Ref. [25| . The weights found 
in this article are small (they are proportional to a) but 
they converge to a constant, as we saw in the paragraph 
above. We use the formula here nonetheless. 

CONJECTURE 5. The critical density of the model 
H5 is given by Eq. (|60[) . 



15. The critical density. The critical density is 
given by Eq. (|60[) . where the cycle weights are given by 
Eq. gHD. To first order, 



y 1 ' i > 1 



j 



3/2 



-C(|) 



(61) 

It turns out that the sums can be explicitly calculated. 
We have 



with 



-1 



2a 



lim Cl 



Cl -E^r/2"[5E( 



(4tt^) 2 l^oo 



3/2 



c(i: 



Interchanging the sums over £ and j , we get 

l . L 



(62) 



(63) 



Cl - E ~m [5 E 

j 1 j t=j+ 



i=j+ 
L-3 



(£ - 7W 2 ~ ^ W 2 
1 v J ' t=i 



E ? -3/2 [2 E 



j > 1 



£3/2 



E li 



£1/2 



(64) 



1 L 

E -3/2 E 
1 ; ^ 1 •> d — 1 T A\ 



1 



3>1 J «=(L-i)++l 

The first line of the latter expression is equal to 



2 E o3/2 E 



1 



1 



j> 1 



^£3/2 £1/2 



= h E 



1 



1 



jl/2£3/2 j3/2£l/i 



(65) 



= 0. 
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We are left with the last line of Eq. (l64l) . Splitting the 
sum over j, we get 

^ L ~ - E -3/2 E 71/2 — E '3/2 E 71/2 



e 1 ' 2 

j=l J £=L-j+l j > L + l J 1=1 

-c (1) - cf ] 



(66) 



After a change of variables, the first term is given by 

L 



C 



(i) 



(67) 



As L — > oo we get a Riemann integral, namely, 



lim C 



(i) 



1 dx 



i.~x ./o x 3 / 2 J (l-x + y) 1 / 2 



dy 



1 1 - VT - ^ 



c 3/2 



dx = -4 + 4 



o - jr) 



(68) 



The last identity follows from an integration by parts. 
The substitution x = sin 2 9 eventually gives 



lim C { P = 2tt - 4. 



The second term in (l66l) is equal to 



(69) 



^tEtT^^E^)- (70) 



I V l - 

L ,rl( 1 +f) 3 



Then 



lim C< 2 > = ' dX fl dy 



-2= =4. (71) 



We have obtained that lirriL Cl = — 2ir, and so we have 
from dBTJ that 



p^ 0) ~c(§r 



/2 



(72) 



16. The critical temperature. We have found 
the constant 6 for the linear change in the critical density, 
i.e. 



(a) (0) 
Pc - Pc . L _o_l/2 



Pc 



= bafi- 



(73) 



We would rather know the constant c for the linear 
change in the critical temperature, i.e. 



-,(a) ^,(0) 



(0) 



= cap 1 / 3 . 



(74) 



The critical line can be seen as a manifold in the space 
(p, P, a); at a — we have the relation 



C(|) 

(4^)3/2 ' 



Recall that 



da dp dp 
dp dp da 



We have 



da ^^ pP-pM yi ^ 1 



3 C(f) 



dp 

df_ 

dp 2 (4tt) 3 /2^5/2 

dp 1 dT 
0a 



bp-^p 



(75) 



(76) 



(77) 



(78) 



T> „ ~c P ^T = -cPp^. (79) 



Inserting these expressions into (|76")) . and using ([75]l. we 
find 



3C(|) 1/3 



(80) 



We can now compute the constant for the change in 
the critical temperature. Using b = 2y / 7r/£(|) from Eq. 



(81) 



72]) and the relation above, we finally obtain 
8tt 



3C(3/2) 4 / 3 



-2.33. 



The negative sign means that interactions discourage 
Bose-Einstein condensation. Our finding is in contra- 
diction with the physicists' consensus c ~ 1.3. 

17. A test OF the method. We now compute the 
free energy of our approximate model, and we compare it 
with the formula for the dilute Bose gas, which is known 
explicitly [19, Ell . 

The free energy of the model H4 is defined by 



/CM 



1 



•log 



P\M (47T/3) 



\ N y [ dxc- H ^. 



(82) 



It is related to the free energy / of the model H$ by 

/(/?, p) = Airap 2 + Sirap^p + f(ft, p). (83) 

The existence of the thermodynamic limit and the equiv- 
alence of ensembles were shown in Ref. (25[. Then / is 
given by the Legendre transform of the corresponding 
pressure, 

f(P, p) = sup [pp -p(P,p)]= pp* -p(P,p*)- (84) 

The pressure for H$ has already been computed, the re- 
sult can be found in Eq. (I5T)1) . To first order, 



p(p,p)=p^(p,p)-J2 



' (4ir) 3 / 2 (p£) 5 / 2 ' 



(85) 
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The equation for p* is p = dp/dp; for p p c , 

4 Q )(R i - \ 11L. (86) 



P = p-H(3,p*)-J2 



' (4tt/^) 3/2 ■ 
Notice that p* = when p ^ p c . We have 



e^* e a' 



' (4tt) 3 / 2 (/3£) 5 / 2 



(87) 



Since p* differs from the maximizer of pp — p^ by 0(a), 
we have 



pp* -p<- \l3,tx*) = / luJ (Ap) + 0(a 2 ). (88) 



f(0) 



We now calculate the last term of (|87[> ■ Using the ex- 
pression |46|) of o^, we get 



E 



2o 



' (47r) 3 / 2 (/3£) 5 /2 /3(4tt/3) 2 

2 7 '3/2f/ _ ? -)3/2 ^l) £3/2 

£ ^ 1 j'=l O 1 

(89) 

The first term of the right side is equal to 47rap(°)(/3, p*) 2 , 
and the second term is equal to 8nap^ p(°\p, p*). We 
have 

(0)(/3 , } fp + 0(a) ifpJPc 

p c + O(O) II P ^ Pc- 



The right side is indeed equal to the free energy of the 
dilute Bose gas up to a correction Q(a 2 ). This formula 
was rigorously established recently [20, [21] . 

The test of the free energy does not prove the validity 
of our method, but it seems fairly robust nonetheless. It 
would detect a mistake in the definition or in the cal- 
culation of the weights. We can also be confident that 
no important term was left behind when the model was 
progressively simplified. It would be interesting to apply 
it to the other methods that served to compute the con- 
stant c. Our understanding would be greatly enhanced if 
we knew which methods pass it and which fail to pass it. 



18. Conclusion. The main goal of this article is to 
propose a method that allows to compute the effects of 
interactions on the critical temperature of the dilute Bose 
gas. We used the concept of "spatial random permuta- 
tions", whose main advantage is to allow the heuristics 
and the intuition of statistical mechanics. Our findings 
contradict those of our colleagues from physics and this 
casts a shadow on the method. We have been careful 
to state our conjectures precisely, which will help resolve 
the discrepancy in the future. We have also performed a 
test (the calculation of the free energy) and the method 
passed it. Besides the main goal, our approach suggests 
several models that may be of interest as effective models, 
even if they do not describe the dilute Bose gas exactly. 



Then 

M^/ ( °»(^) + ( 4 7 (p ; 2pc) :| p ; pc (3D 

y-^irap* it p > p c 

Together with ([55]) . we get 

/CM = / (0) (/3,p) +47ra[2p 2 - (p - p c )\\. (92) 
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